On One Model of the Geometrical Quintessence 
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C/3 I A cosmological model with modification of the Einstein-Hilbert action by the 

■ correction f(R) oc (3K n is considered. Such way of the description of the "geo- 

metrical" dark energy has been introduced repeatedly and the coefficients of the 
model were chosen to be fit against some observational data. In this paper the 

T — I ' unambiguous choice of parameters n and (3 is proposed from the follow reasons: the 

exponent n close to 1.2953 follows from the request for the evolution of the Universe 
\ after recombination to be close to the evolution of the flat FRW model with cold 

dark matter and the reasonable age of the Universe defines the magnitude of the 
coefficient j3. Such a model corresponding to the evolution of the Universe with the 
dynamical A-term describes well enough the observational data. 

1 Introduction 

£3 ■ 

The discovery of accelerated expansion of the Universe [B El El has stimulated the quest 
>• . for mechanism of the modern inflation. The most famous theoretical model of dark energy 
(DE) is the cosmological constant A. The corresponding FRW solution for flat Universe 
with the present densities ratio for cold matter and dark energy (Q m /Q\o ~ 0.3/0.7) 
describes satisfactorily the evolution of the Universe at low redshifts jSHHEJ- However, 
the nature of the constant A- term has been remaining to be inexplicable during many 
years. 

It is well-known, the application of the constant A- term for modeling of the early 
Universe was initially confronted with principal difficulties. Solving the problems of the 
very early Universe, this term has to be reducing by several order of magnitude during 
following evolution of the Universe. This problem was solved by rejection of A- term, and 
corresponding inflationary behavior was determined by models of " effective A- term"- 
quasi-classical scalar fields in the one way. The progress of these models is well-known. 
Another way to describe the inflationary behavior is to take into account the polarization 
of vacuum of quantum fields in the early Universe. Taking into account of the effects 
of the polarization leads to the appearance of the terms non-linear on curvature in the 
Einstein-Hilbert action. In such models the inflation appears self-consistently. Let us 
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note two issues 

- a correction to the Einstein -Hilbert action with the arbitrary function of scalar curva- 
ture R is equivalent mathematically to the introduction of scalar field into the classical 
Friedman cosmology jfH 13 IE) ; 

- the terms of the form R n were investigated in the early works done on the problem of 
singularity before obtaining of exact corrections to Einstein action following from the one 
loop approximation j§J . The part of such solutions approaches asymptotically to Friedman 
solutions with A-term, however physical results of the solutions have not been explained 
at that time. 

For the purpose to explain the accelerated expansion of the Universe today, it is natu- 
rally to use the experience cumulative in investigations of early Universe. Thus, one of the 
tendencies is concerned with the hypothesis of existence some scalar fields that determine 
the density of dark energy (DE) (e.g. [TU1 ITT1 IT2*]). The another tendency models an 
effective quasi-hydrodynamical tensor of momentum-energy describing the observational 
data ^3]. And third tendency is to generalize the Einstein-Hilbert equation by inclusion 
of curvature invariants ED E] analogously to the more early works mentioned above. 
The last approach one can consider to be either an independent approach to the describing 
of DE or an analogue of inclusion of scalar fields (in the case f(R)) as stated above. 

In the works on the higher order gravity theories (HOGT) the models with power 
corrections were investigated, however they have never been fitted to whole set of the 
observational data. 

In this paper, the model with correction f(R) oc R n with n > is considered in 
detail for the purpose to correlate it with the observational data. In the other words, 
we would like to obtain the model that does not conflict with the scenario of the large 
scale structure formation (in past) and describes satisfactorily the Universe undergoing 
an accelerated expansion at present. Therefore, at the minimum, in the framework of 
the f(R) -theories, we will obtain solutions remind ACMD model describing well enough 
by set of the observational data Ej. However, as it has been mentioned by various 
authors (see [5] and references therein) the observational data indicates the models of the 
dynamical DE. Hence, our second aim is to search out such dynamical solutions in the 
framework of the HOGT and to find out whether these solutions are preferable. 

This paper is built as follows. In the section II the basic equations of the HOGT 
are presented. In the section III for the corrections of the form f(R) oc (3R n we find the 
exponent n\ = 1.2953 which allows the generalized Einstein equation for the scale factor a 
to have a particular solution corresponding to the flat FRW solution for cold matter. We 
show that instability of the solutions that are close to this particular solution at z » 1 
may lead to the accelerated behavior of the model at present and the following asymptotic 
approach of the solution to the solution with the constant A- term. In the section IV we 
discuss our results and compare them to observations. In the model there are only two 
free parameters of the model - the coefficient (3 and a slight deviation of the parameter 
n from n\ = 1.2953 mentioned above. Fixed from one set of the observational data they 
allows to obtain the rest of the set of the observational data. 
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2 Basic equations 

We will work in the matter frame with a spatially flat FRW metric 

ds 2 = dr 2 - a 2 (r)(dx 2 + dy 2 + dz 2 ), 



with a(r) being a scale factor. The non-dimensional components of Ricci scalar and Ricci 
tensor are, 



R=-6 



+ 



Rq - -3 



1 da 
H dr 



(2) 



Here the present Hubble parameter Hq has been introduced to turn to non-dimensional 
values. The equations that can be derived by variation of the action 



S=^fj d 4 xV^ (R + f(R)) + J d 4 x^ L r 



(3) 



are well-known for a long time [TJj- It is useful to reduce its — component , that is the 
third order equation in a, to the second order equation by introduction of a new function 



y(a) = (da) . 

Thus, Eqs.(|2|) are determined as 

R = -3y'/a 2 , R% = 3(y - ay')/2, y' = dy/da, 



(4) 



(5) 



and the observable parameters - the Hubble parameter h and the deceleration parameter 
q - are determined by formulae 



h(a) 



aa 



1 



dHy) \ 
d\n(a 2 ) J ' 



(6) 



Further, the — component of Einstein equation can be represented in the form 



y + 



df 



a dy 



d f df 



dR 



2 da 



6 



da \ dR 



a. 



(7) 



The solution for the classical Friedman model with cold matter is defined by f(R) = 
whereas for the model with the constant A-term it is defined by f(R) = const. Taken 
together they describe ACDM model, so, according to Eqs.©,© we have 



y = fl m a + fl\ a 



n, 



Pm/Pc, Qa = Pa/Pc 



where p m , p\, and p c = 3Hq/8ttG are the cold matter density, the dark energy density, 
and the critical density of FRW Universe at present, respectively. For 4- interval Q 
fl m + Q\ — 1. The present-day observational values are close to Q m = 0.3, = 0.7. 



3 



Notice that the important property of Eq. (|7jl is a possible asymptotic exit of solutions 
to regime of an exponential expansion. Indeed, according to definition JHJ) this solution 
is y = Qa 4 . In this case, the scalar curvature approaches to constant Roc —>■ —12Q. In 
the case of a negligible effect of matter this solution is a particular solution of Eq.fJZj) at 
execution of the following condition 

3 Model with Power Corrections in Einstein -Hilbert 
Action 

In this section we will investigate corrections to the Einstein-Hilbert action of the form 

f(R) = -aR n . (10) 
In such Eq. (0J) can be presented in the form 



P 



n{n - l)y"y + [ ^-^{y') 2 + n(4 - 3n)^ 

— (X 



(y - il m a) 



a 4 " 3 " 



/3=(-3)"- 1 a. (11) 



A general approach to investigation of the last equation is given in paper 118] . 

There are two parameters in the model, f3 and n, determined by observational data. In 
a number of papers these parameters were chosen according to different requirements ^Hl 
ITBj . In the present paper we will choose parameter n from requirement of closeness of 
evolution of the our model to the classical solution for the flat FRW Universe with cold 
matter in the past. This fact allows this scenario to be close to the scenario of the large 
scale structure formation. This requirement can be realized on condition the classical 
Friedman solution 

y = Q m a (12) 

presents a particular solution of Eq. lfTTj) . It is easy to see from Eq. ffTTj) . the last condition 
is equivalent to the choice of n to be satisfying the equation 

n - 1 = 2n(4 - 3n) with roots m = 1.295, n 2 = -0.129. (13) 

The first of the roots leads to the type of models of papers (13 EH , while the second 
root corresponds the models with correction of the form oc /x/i?'™ 2 ' investigated in the 
paper [20] • As it will be shown below, such a choice of n approaches the model to the set 
of the observational data in the best way. 
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3.1 Behavior of the dust solution in the /(i?)-model in past 

After recombination, the evolution of the Universe has to be described by the Friedman 
model with the cold matter. The dynamics of expansion is determined by stability of the 
dust solution in the model (jllj) . If the solution is stable, then the model evolves in the way 
very close to classical one. However, the observational data at z < 1 does not correspond 
to such a scenario, i.e. we are interested in dust-like solutions which are not stable in 
the model (|TT|) but the perturbations do not grow catastrophically fast. Otherwise, the 
model does not provides a sufficiently long period with g~0.5at2^>l required for the 
large structure formation. 

For investigation of behavior of the solution (j!2|) we will search a perturbed solution 
in the form 

y = fia(l + V0, ^<1- (14) 
In this linear approximation and close to the recombination time (a <C 1) Eq. fTTf yields, 

n a 2 ip" + a ip\2n - 0.5) = 0, (15) 

with the damping solution ip = Ci + C^a 1 -^ -1 ) . Hence the solution ()12|) for /(i?)-theory 
(fTUj) with n = rii satisfying Eq.ffT^j) asymptotically approaches to to the flat FRW solution, 
i.e. it is stable and does not satisfies the requirement stated above 1 

As the next step on the path of the choice of n, we will look for solutions of the f{R) 
-theories (|l(Jj) with n which is a little different from ni, 

n = rii + 5n, 5n <C 1. (16) 

It is efficient to rewrite the condition (|13|) in the form 

5(n- 1) = 2n(4-3n), 5 = 1 + e, e < 1. (17) 

We will look for the perturbed dust-like solution in the form (|14j) near to recombination. 
The Eq. (fTT|) for such a case yields 

n a 2 ip" + 2a ip' (n - 0.5 + 5 /A) + (5 - 1) = 0. (18) 

Then after the change of variable £ = ln(a/a*) the Eq. (|T5|) yields 

nip +{n- 1/2 + e/2)ip + e = 0, m =d/d£. (19) 

J It is interesting to note that the mentioned exponents (|1 H|) are obtained in the recent paper ^1 fc> r 
the equation equivalent to Ea. (|ll|l . In f(R) theory, the given exponents allow obtaining of the solutions 
for cold dark matter coinciding with ones in the classical Friedman model of the universe. Also, it has 
been shown in the paper ^] that these solutions are stable in the framework of HOGT. Let us notice 
we slightly change exponent n in the preset work to obtain weakly unstable solutions adequate to the 
observational data. Authors are grateful to authors of paper JU] kindly attracting our attention to their 
results. 
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Figure 1: The evolution of the variable y (a left panel) and the Hubble parameter h (a right panel) with 
the scale factor a are represented by the solid line. The ACDM solution and the cold matter solution are 
presented by dashed and dotted line, respectively. Both graphs correspond to the case O m = 0.29. 

The last equation has a solution 

M^ + T^ 1 ' 1 ©^^ 1 ' (20) 

The analysis of (|2(J|) has shown that the requirement stated above is realized only for 
e < 0. 

In this case, the modification of the exponent (fTTjl is determined by a small positive 
correction 

5n~ -2en x {n x -\)j{Yln x - 1). (21) 

The numerical analysis has shown that behavior of solutions is sensitive to small changes 
of e at a — * 1. The last fact together with a choice of parameter j3 allows us to obtain a 
good enough correspondence with the observational data. 

As an illustration we give the results for the set of parameters n = 1.296, (3 = 0.467 
fixed according to Q m = 0.29, h = 68 (the best fit to the CMB+SNe data presented 
in the paper Joj. In a left panel of Fig. Q a solid line represents the evolution of the 
variable y with a scale factor a. At the beginning, it coincides with the evolution of the 
dust model which is represented by a dotted line but further it deviates to the ACMD 
model represented by a dashed line. In a right panel, one can see evolution of the Hubble 
parameter with the scale factor a for the mentioned three models. 

3.2 The behavior of the solution in future 

The further expansion of the Universe at a » 1 according to the Eq. (|11|) leads to the 
negligible effect of cold matter on the solution behavior. In this case, De Sitter solution 



6 



30 

a 



60 



30 

a 



Figure 2: The evolution of the Hubble parameter h and the deceleration parameter q with a scale factor 
a. The thick solid lines represent graphs related to the /(i?)-model while the dashed lines represent the 
evolution of the parameters of ACDM model. In the left panel the asymptotic solution for f(R) model 
is shown by the dash-doted line. 

y = Qa A is an asymptotic solution of Eq. fjllj) . This solution corresponds to the constant 
Hubble parameter 

/i(a^oo) = Vti (22) 

with Q defined from equality 

3/3(2ra- l)(n- l)/2 = (4fi) (1 - n) /4 (23) 

The inflationary solution is stable in the process of evolution of the model. To show 
it, we shall look for the solution with perturbation in the form 
y = fia 4 (l + $), $ << 1. This ansatz yields 

niaV + [(ni - l)/2 + 3n?]a$' + 6(2n? - 3m + 1)$ = (24) 

The change of variable a to variable £ (see Eq.fl^UJ)) yields 

<£ + A$ + fi$ = 0, ' = d/df . (25) 

where coefficients are A = [3m - (1/m + l)/2] = 2.99, 5 = Q(2n x - 3 + 1/m) = 2.18. 
This equation for perturbations have a damped solution indicating De Sitter solution 
to be stable. The numerical analysis has shown that De Sitter solution is an attractive 
solution. 

The evolution of the Hubble parameter is represented by the thick solid line in the left 
panel of Fig. El for the case fl m = 0.29. In contrast to the Hubble parameter of ACDM 
model (the dashed line) monotonically decreasing down to constant it reaches a 



minimum h m ~ 0.978 at a ~ 1.15 and after that increases up to the asymptotic solution 
(122}) which is represented in the figure by the dash-doted line. It is interesting to note that 
the formula for dimensionless parameter Hubble h(z) = h(a) obtained from observational 
data in the paper jH] allows it's extrapolation in future (a > 1). At the parameters men- 
tioned in this paper, the formula for h(a) also predicts minimum of the Hubble parameter 
at a ~ 1.45 which is equal to 0.951. This best fit of the paper jS] is shown in Fig|2]by the 
thin solid line. The deceleration parameter q represented in the right panel of Fig. |21 also 
passes a minimum and approaches to —1 with the growth of a. Therefore, we live in a 
transitional epoch between the classical Friedman cosmology and the De Sitter cosmology. 

4 Results 

Hereafter we will present the comparison of results of our model and the observational 
data. 

As a set of observational data, the analysis of SNe+CMB data done in the paper [5J 
have been used. In that paper authors have reconstructed the resent history of the 
Universe on the base of SNe and CMB data in the model-independent way, only modelling 
DE by the hydrodynamical equation of state 

p = wp, w = (2q-l)/(3-n m /h 2 ). (26) 

The cited paper presents two conceptions of the analysis of observational data: the first 
of them is the best fit to the data which uses only the hydrodynamical describing of DE 
and does not impose restrictions on the values of fl m and ho, while the second conception 
follows the priority of the concordance ACDM model, so authors of [3] put fl m = 0.27±0.04 
and h = 0.71 ± 0.06. 

We will give the comparison of our results with both of them. Also we notice that 
the analogy of the "hydrodynamical" DE (|27)|) is not so proper to the higher order grav- 
ity theories, hence one can expect the comparison over the values h and q to be more 
informative than over the value w. 

It have been found in the paper [5] the best fit values are: Q m = 0.385, h = 60. In the 
model DE evolves in time strongly enough. In Fig. El one can see results (the evolution of 
the deceleration parameter and DE equation of state parameter with redshift) for analysis 
of paper |S]: the best fit is represented by the thin solid line, the la and 2a confidence 
levels are represented by the light and dark grey contours, respectively, and ACDM is 
represented by the dashed line. Also for given fl m and h the results of the /(-R)-model 
with n = 1.2955 (5n = 0.0002) and (3 = 0.273 for the "geometrical equation of state" 
parameter w = (y — ay')/3(y — Q m a) and deceleration parameter q are presented by the 
thick solid lines. In given /(i?)-model the age of the Universe is 14.9 Hyr, the deceleration 
parameter is go = —0.91 at present and the transition to acceleration occurs at z = 0.38. 
Similarly to results of pj, the wde < 1 at lower redshifts (wdeo = —1.53), however , the 
evolution of equation of state of "geometrical DE" is more weak contrary to the results 
of 0. 
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Figure 3: The comparison of results of the f(R) oc f3R n model to the analysis of the observational data 
quoted from the paper [5] by approbation of the authors. The evolution of the deceleration parameter 
with redshift is shown in a right panel and the variation of the parameter w of equation of state of DE 
is shown in a left panel. The best fit of SNe+CMB data with f2 m = 0.385, ho = 0.60 is represented by 
the thin solid line, the la and 2a confidence levels are represented by the light and dark grey contours, 
respectively, and ACDM is represented by the dashed line. For given fl m and h the /(i?)-model is defined 
by parameters n — 1.2955 and (3 — 0.273 and represented here by the thick solid line. 




Figure 4: The comparison of results of f(R) oc R n model with n = 1.296 and (3 = 0.467 (the thick 
solid line) to results of analysis of SNe+CMB data with f2 m = 0.29 done in 0. The evolution of the 
deceleration parameter with redshift is shown in a right panel and variation of equation of state of DE is 
shown in a left panel. The best fit of SNe+CMB data in such case is represented by the thin solid line, 
the la and 2a confidence levels are represented by the light and dark grey contours, respectively, and 
ACDM is represented by the dashed line. 
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However, if strong priors have been imposed on Q m and ho (i.e. the ACDM model 
priors: Q m = 0.27 ± 0.04 and h = 0.71 ± 0.06), the evolution of DE is extremely weak 
and in good agreement with the ACDM model. The best fit in the case is Q m = 0.29. 
One can see good enough coincidence of our model and their analysis for parameters of 
the model n = 1.296 (5n = 0.001) and (3 = 0.467 in figure |U The deceleration parameter 
q = —0.683 at present and the deceleration was changed by the acceleration at z — 0.51, 
(go = —0.63 ± 0.12 and z = 0.57 ± 0.07 in the paper [3]). The age of the Universe in this 
case is 13.6 Hyr. 

In Fig. El the comparison of the evolution of the Hubble parameter for /(i?)-model to 
ACDM model and the best fit of the observational data analysis is presented. One can 
see the results for the cases with Q m = 0.385 and Q m = 0.29 in a left and a right panel, 
respectively. 

Thus, the f{R) oc /3R n -mode\ with parameters (3 and n chosen according to the princi- 
ples mentioned in Introduction describes the evolution of the Universe quite corresponding 
to the SNe+CMB data. 




Figure 5: The evolution of the Hubble parameter with redshift. The /(i?)-model and the ACDM model 
are shown by the thick solid and the dashed lines, while the best fit of the observational data according 
to p] is represented by the thin solid line. The case for fl m — 0.385, ho — 0.6 is shown in the left panel, 
and the case f2 m = 0.29, ho — 0.68 is shown in the right panel. 
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